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^ . Abstract. Here we report on measurements of the absolute absorption and dispersion 

£>~V properties of an isotopically pure 87 Rb vapour for magnetic fields up to and including 

■ 0.6 T. We discuss the various regimes that arise when the hyperfine and Zeeman 

interactions have different magnitudes, and show that we enter the hyperfine Paschen- 
Back regime for fields greater than 0.33 T on the Rb D2 line. The experiment uses a 
compact 1 mm 3 microfabricated vapour cell that makes it easy to maintain a uniform 
I and large magnetic field with a small and inexpensive magnet. We find excellent 

C ■ agreement between the experimental results and numerical calculations of the weak 

probe susceptibility where the line positions and strengths are calculated by matrix 
diagonalization. 
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1. Introduction 

The interaction of light with atomic ensembles continues to be a topic of active research. 
Achieving strong coherent matter-light interfaces allows us to investigate applications 
within quantum information processing, such as the quantum internet pQ and quantum 
memories [2]. The addition of an external magnetic field for such an ensemble further 
increases the possibilities for these interfaces to be realized. Resonant and off-resonant 
linear and nonlinear magneto-optical effects in multi-particle ensembles have allowed 
the observation of quantum teleportation between light and matter [3] and control of 
atomic Zeeman populations in long-lived quantum memories [I]. 

As first observed by Zeeman in 1896 [5], the energy levels and transition 
probabilities of an atomic ensemble are extremely sensitive to an external magnetic 
field. Therefore a theoretical model for the absorption and dispersive properties of 
such ensembles has found utility in many applications, for example, Faraday dichroic 
beam splitter for Raman light [6]; Gigahertz-bandwidth atomic probes [Zj; Hanle-type 
coherent population trapping [8]; off- resonance laser frequency stabilization [9]; realizing 
narrowband atomic filters [10J; cooperative effects in an atomic nanolayer imaging 
microwave fields in vapour cells p2] and achieving a compact optical isolator [13] . 

The absorption and dispersion of an atomic ensemble in an external magnetic field 
can be calculated from the atomic susceptibility. In the absence of field, absolute 
Doppler-broadened absorption [Til [15] and dispersion [16] in the low density regime, 
and dipole-dipole interactions [17] in the binary-collision regime have been tested. In 
the presence of field, Stokes parameters for fields up to 0.08 T [18] have also been 
investigated. The motivation for this work is to test the model for absolute susceptibility 
in a 87 Rb vapour on the D2 line. To the best of our knowledge this is the first 
study highlighting excellent agreement between experimental results and numerical 
calculations for absolute absorption and dispersion for fields up to and including 0.6 T. 
It should be noted that at such fields the nuclear spin, /, and total electronic angular 
momentum, J, are decoupled for both ground and excited states and the total angular 
momentum, F, is no longer a good quantum number; this is known as the hyperfine 
Paschen-Back regime [19]. 

The structure of the rest of the paper is as follows: in Section [2] we discuss the 
quadratic magnetic term for fields encountered in the laboratory for low- and high-lying 
n states and also the different regimes available for the linear magnetic term as a function 
of magnetic field. In Section [3] we describe the isotopically pure 87 Rb microfabricated 
vapour cell and permanent neodymium magnet used in the investigation. In SectionHJwe 
measure the evolution of the absolute optical depths as a function of field and detuning. 
To highlight the decoupling we also show the transition-strength dependence for the 
outermost weakly allowed transitions. In Section [5] we compare theory and experiment 
for the mediums absolute absorption and dispersive properties in the hyperfine Paschen- 
Back regime, highlighting the excellent agreement. Finally, in Section [BJ we draw our 
conclusions. 
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2. Theoretical Considerations 

2.1. Atomic Hamiltonian 

The atomic Hamiltonian can be written as 

H = Ho + H is + H h{s + H z , (1) 

where H Q is the coarse atomic structure; H fs and H h f s describes the fine and hyperfine 
interactions and Hz] represents the atomic interaction with an external magnetic field. 
The zero detuning frequencies in the absence of hyperfine splitting for 87 Rb and 85 Rb 
on the D x (5 2 S V2 5 2 P 1/2 ) and D 2 (5 2 S 1/2 5 2 P 3/2 ) lines are 377.11 THz [20J and 
384.23 THz [21], respectively. The splittings associated with the hyperfine interaction 
Hamiltonian around the zero-detuning energies can be calculated by use of the following 
expression 

AP A hfs B his \K{K + 1) - 21(1 + I) J {J + 1) 

AEhfs = — K + 1 J(2/-1)J(2J-1) ' (2) 

where A Ms is the magnetic dipole constant, B h f s is the electric quadrupole constant and 
K = F(F + 1) — 1(1 + 1) — J (J + 1); see equation 9.60 in [22]- The numerical values to 
this expression can be found in, for example, table 1 of [14]. For an external magnetic 
field the magnetic interaction Hamiltonian has a linear and quadratic term in 3 [23J. 
For typical laboratory magnetic fields, quadratic shifts for states with a low principal 
quantum number, n, are extremely difficult to observe. Therefore the quadratic term 
in B can usually be ignored; however, with high n states the dependence is very much 
evident 0- The magnetic interaction Hamiltonian for an external field has the form 

#z = -(Ai, + Mz + A*s)-B, (3) 

where fi L and fi s are the magnetic moments of the nucleus, the orbital motion due 
to the electron, and the spin of the electron, respectively. We can ignore the contribution 
due to the magnetic moment of the nucleus, /x j5 because the Bohr magneton is three 
orders of magnitude larger than the nuclear magneton [23]. The different spacings 
AEz associated with the magnetic interaction Hamiltonian are discussed in section 12.21 
In this work we investigate two cases: the hyperfine linear Zeeman regime (HLZ), 
where the magnetic interaction is treated as a perturbation to the hyperfine interaction, 
AEz < AEhf s ; and the hyperfine Paschen-Back regime (HPB) where the magnetic 
interaction is larger than the hyperfine interaction, yet smaller than the fine interaction, 
Ai^hfs < AEz- For even larger fields the magnetic interaction dominates both the fine 
and hyperfine interaction; this is known as the fine Paschen-Back regime (FPB). 

2.2. Various magnetic regimes 

Figure [1] shows the energy level evolution as a function of magnetic field for 87 Rb 
(/ = 3/2) on the D 2 line. In the weak-field regime the nuclear spin, /, and total 

\ The study of diamagnetic shifts in Rydberg atoms has been studied for many years. For magnetic 
fields of about 0.6 T the shifts become evident for states above n ~ 36 [24] . 
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Figure 1. Diagram showing the shift in the energy levels as a function of magnetic 
field for the 5 2 Si/ 2 and 5 2 P 3 / 2 terms in 87 Rb. In the (weak- field) hyperfme linear 
Zeeman (HLZ) regime, the energy levels can be described by the F and mp quantum 
numbers (shown in the left plot). At intermediate fields there are no good quantum 
numbers to label all of the energy levels. At large fields (0.6 T) mj and mj become 
good quantum numbers due to the nuclei and electronic spins decoupling (shown in 
the right plot); this case is referred to as the hyperfme Paschen-Back (HPB) regime. 



electronic angular momentum, J, couple to give the total angular momentum, F, which 
has 2 J+l values for J < 1. For the 5 2 S"i/2 term of 87 Rb, F can be 1 or 2 with a hyperfme 
splitting of 6.8 GHz, whereas for the 5 2 P^/2 term, F can be 0, 1, 2 or 3 with hyperfme 
splittings between 70 and 270 MHz. For weak magnetic fields each hyperfme level, F, is 
split into 2F+1 levels (mp) symmetrically about the zero field level, AEz = QfTtif^bB] 
this is known as the HLZ regime. In this regime the \F,mp) basis best describes the 
interaction. For large magnetic fields the total angular momentum decouples into I and 
J; this is known as the HPB regime. The effect introduces 21 + 1 levels (mi) with each 
mj value, for the 5 2 Si/2 term, mj is equal to 1/2 or —1/2 and for the ^> 2 Pz/2 term, mj 
can be 3/2, 1/2, —1/2 or —3/2. In this regime the \mj,mj) basis best describes the 
interaction. The spacings of the levels are proportional to the values of rrtj and mj, 
AEz = (gjirij + Qittl^ubB. A detailed theoretical discussion of this regime can be 
found in [23]. For intermediate magnetic fields all symmetry is lost and there is no good 
basis to describe the interaction. For even larger fields, typically > 218 T, J decouples 
into the total orbital angular momentum, L, and total spin angular momentum, S; 
this is known as the FPB regime. This effect introduces 2L + 1 levels (mj) for each 
of the 2 orientations of mj for a single electron. These are known as Lorentz triplets 
when L — 1. The spacings of the levels are proportional to the values of mx and ms, 
AE Z = (gLm L + g s ms)fi B B. 

Table [1] shows the typical calculated magnetic fields to gain access to the HLZ, 
HPB and FPB regimes of the linear magnetic interaction term of 87 Rb and 85 Rb. In 
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Table 1. The magnetic fields required to gain access to the linear and quadratic 
terms of the magnetic interaction Hamiltonian for Rb. For a field of 0.6 T and n = 5 
state only the hyperfine linear Zeeman (HLZ) and hyperfine Paschen-Back (HPB) 
regimes are accessible in this work. The magnetic dipole constants, A^s, and electric 
quadrupole constants, £?hfs, were obtained from [21 . The fine Paschen-Back (FPB) 
regime is also shown. The n w 36 state would need to be investigated to see any 
quadratic effects for the field under investigation, or a field of w 1.6 x 10 3 T would be 
required to investigate the n = 5 state. 







Linear 


Quadratic 


Isotope 


Term 


HLZ / HPB 


FPB 






A hk /h (MHz 


) £? hfs /h(MHz) 


B(T) 


B(T) 


B(T) 




5 2 S 1/2 


3417.34 




0.33 






87 Rb 


5 2 Pi/ 2 


407.24 




0.12 


218 






5 2 P 3/2 


84.72 


12.50 


0.01 


« 1.6 x 10 3 




5 2 S 1/2 


1011.91 




0.13 




85 Rb 


5 2 Pi/ 2 


120.32 




0.05 


218 






5 2 P 3 /2 


25.00 


25.79 


0.005 






n ~ 36 






0.6 



addition, the fields and principal quantum number required to see any quadratic effects 
are noted. To calculate the fields we equate the hyperfine energy splitting to the sum of 
the Zeeman shift of the lowest rriF state of the upper F value and the highest rriF state 
of the lower F value. The HLZ and HPB regimes are accessed for fields much smaller 
and bigger than the calculated values, respectively. A similar procedure is adapted to 
calculate the fields required to access the FPB regime. For the low-lying n = 5 state, 
fields of ~ 1.6 x 10 3 T are required to observe the quadratic term of the magnetic 
interaction Hamiltonian. For 0.6 T we would need to gain access to the n « 36 state 
to measure any such interaction. 

2. 3. Matrix representation of the Hamiltonian 

A detailed description of the model used to calculate atomic susceptibility incorporating 
the energy levels and transition probabilities of the Rb ensemble can be found in [18J. 
In summary the lineshape around resonance is given by a convolution of the Lorentzian 
(accounting for natural, self-broadening and buffer gases) and a Gaussian distribution 
incorporating the Doppler shift due to thermal motion. The total susceptibility, \i 
is then calculated by summing over the electric-dipole-allowed transitions. From the 
atomic susceptibility we are able to model the absorptive and dispersive properties of the 
medium. The refractive index and absorption coefficients can be calculated by use of the 
real, n = l + 9ft(x)/2, and imaginary, a = k^s(x), parts of the susceptibility, respectively, 
where k is the wavevector. To obtain numerical values for the susceptibility a matrix 
representation of the Hamiltonian for the fine and hyperfine interactions and atomic 
interaction with an external magnetic field are calculated in the completely uncoupled 
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\mj, mi, m s) basis. The frequency detunings and transition strengths are calculated from 
a numerical diagonalisation of the matrix. 

3. Experimental Details 

The experimental procedure used to investigate the HPB regime is described in [13] . In 
summary, we use a linear polarized weak-beam |25j with a power of 10 nW and 1/e 2 
radius of 80 /im propagating along the z-direction. The laser traverses a microfabricated 
cell of length L = 1 mm. The frequency axis is calibrated by use of the method described 
in [18]; for detunings over 60 GHz, two scans were stitched together in this experiment. 
Transmission spectra are measured using a single calibrated photodiode. Dispersion 
spectra use a balanced polarimeter to measure the light intensities of the horizontal, 
I x , and vertical, I y , channels after a polarization beam splitter. A half- wave plate is 
set such that in the absence of rotation both channels of light are equal [26]. For this 
publication we describe in detail the two main experimental components: the permanent 
neodymium magnet and 87 Rb microfabricated vapour cell. 

3.1. Permanent neodymium magnet 

To gain access to magnetic fields to investigate the HPB regime, an axial magnetized 
annular permanent neodymium magnet with a circular bore was chosen. The direction of 
the beam in this experiment is parallel to the ^-component of the magnetic field. There 
is an analytic solution for the axial field of a uniformly magnetized annular magnet, 
which is [27] 



where 2t is the length, d = 2r is the inner diameter of the magnet, D = 2R is the 
outer diameter and B r is the remanence of the magnetic material. Figure |2] shows the 
measured solid (blue) circles and theoretical solid (black) line comparison for the axial 
variation of the ^-component of magnetic field. The error on the measured field and 
position is less than the size of the data point. Values of 2t = (6.18 ± 0.12) mm, d = 
(7.98 ± 0.10) mm and D = (25.0 ± 0.6) mm were extracted from a Marquardt-Levenberg 
fit [28] by allowing the dimensions of the magnet to be free parameters. The parameters 
are consistent with the physical dimensions of the magnet, with the deviation between 
theory and experiment at the 0.4 mT rms level. The excellent agreement validates the 
assumption of a uniform magnetization for the magnet. The best-fit remanence for this 
magnet is B r = (1.42 ± 0.07) T. 
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Figure 2. Dimensions of a permanent neodymium magnet and axial variation of the 
z-component of magnetic field. The length, 2t, inner diameter, d, and outer diameter, 
D, describe the dimensions of the magnet, and the remanence, B r , characterizes the 
strength of the material. The measured solid (blue) circles are achieved by use of a Hall 
probe, and the theoretical solid (black) line is obtained from equation 21 Below the 
main graph is a plot of the residuals (solid red circles), that show excellent agreement 
between theory and experiment, with an rms deviation of 0.4 mT. From a Marquardt- 
Levenberg fit, the parameters were found to be: 2i = (6.18 ± 0.12) mm, d = (7.98 ± 
0.10) mm, D = (25.0 ± 0.6) mm and B r = (1.42 ± 0.07) T. 

3.2. 87 Rb microfabricated vapour cell 

A 1 x 1 x 1 mm 3 heated isotopically pure 87 Rb microfabricated vapour cell |29j was 
used that contained buffer gases including hydrogen (H2) and methane (CH4). To obtain 
the optical depths required to compare theoretical and measured absorption spectra, the 
cell was placed in an oven allowing the number density to vary over several orders of 
magnitude. Figure [3] shows a plot of the transmission through the heated vapour cell 
on the Rb D 2 line. The solid (red) and dashed (black) lines show the measured and 
theoretical transmission, respectively, for several temperatures and number densities. 
Based on the discussion in [17] we have previously accounted for a number- density- 
dependent increase in the Lorentzian width due to dipole-dipole interactions; however, 
we must now also include the broadening and shift due to the buffer gases. The collisional 
broadening and shift of the Rb D x and D 2 lines by rare gases have previously been 
measured [30]. For the Rb D 2 line and buffer gases H 2 and CH 4 , the broadening, 
Fbuffer/2vr, equals 26.4 MHz Torr -1 , 26.2 MHz Torr -1 and the line shift, A buffer /27r, 
equals —3.8 MHz Torr -1 , —7.0 MHz Torr -1 , respectively. From the analysis of figure [3] 
we learn that the cell contains 99% 87 Rb, 1% 85 Rb, and a total pressure of ~ Torr of H 2 
and CH 4 . 

Having characterized the magnet and the microfabricated cell we could then 
perform spectroscopy in the HPB regime. 
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Figure 3. Transmission plots for comparison between experiment and theory for the 
Rb D2 line, through a 1 mm vapour cell (99% 87 Rb, 1% 85 Rb) as a function of linear 
detuning, A/27T, for three different temperatures. The solid (red) and dashed (black) 
lines show measured and expected transmission, respectively. The temperatures of 
the vapours were extracted from a Marquardt-Levenberg fit [28] and were found to be 
(60.4 ± 0.2) °C (top), (90.1 ± 0.5) °C (middle) and (127.4 ± 0.8) °C (bottom). The 
broadening and line shift due to the buffer gases were found to be rb u ffer/27r = (23.7 ± 
1.2) MHz and Ab u ffer/27r = (-7.9 ± 2.1) MHz, for all three temperatures, respectively. 
Below the main figure is a plot of the residuals between experiment and theory for the 
middle measurement. There is excellent agreement between the data and model, with 
an rms deviation of 0.5%. There is, however, a small number of glitches due to the 
linearisation of the laser scan being inadequate. 



4. Absorption in the hyperfine Paschen-Back regime 

Figure @] shows the measured absolute optical depths as a function of detuning and 
magnetic field on the D2 line in a 87 Rb vapour. The model for the susceptibility 
yields solid (grey) theoretical lines corresponding to the evolution of the transition 
frequencies as a function of magnetic field. In the zero-field regime, the solid (blue) 
measured absolute optical depths are shown at a temperature of (116 ± 1) °C, the two 
large features are the Doppler-broadened 87 Rb transitions from the F = 1 and F = 2 
states. The other two very small features arise from the 1% of 85 Rb in the cell. In the 
intermediate regime, the solid (olive) measured absolute optical depths are shown at 
a temperature of (116 ± 1) °C and for a field of (0.180 ± 0.001) T. There are many 
spectral features; assigning quantum numbers to the transitions is difficult because there 
is no suitable basis set: the excited terms are completely decoupled, whereas the ground 
terms are only partially uncoupled. 

In the HPB regime, the solid (red) measured absolute optical depths are shown 
at a temperature of (116 ± 1) °C and for a field of (0.618 ± 0.002) T. The spectrum 
remains very rich in structure; however, the |mj, rrij) basis best describes the interaction, 
with the expected 16 strong transitions at this field being clearly visible. Note that as 
the Zeeman shift exceeds the hyperfine interaction, the spectrum becomes symmetric 
with respect to detuning. The other weaker transitions arise as a consequence of 
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Figure 4. Experimentally measured absolute optical depths for the Rb D2 line, 
through a vapour cell (99% 87 Rb, 1% 85 Rb) of length, L = 1 mm, as a function 
of linear detuning, A/2ir, at three different magnetic field values. All three spectra 
were measured at a temperature of (116 ± 1) °C. The solid (blue) measured spectrum 
was taken in the absence of magnetic field. The solid (olive) measured spectrum was 
taken at a field of (0.180 ± 0.001) T in the intermediate regime. The solid (red) 
measured spectrum was measured at a field of (0.618 ± 0.002) T in the hyperfine 
Paschen-Back (HPB) regime. The solid (grey) theoretical lines show the transition 
frequencies as a function of magnetic field. Also plotted is the solid (yellow) theoretical 
transition strength of the outermost weak transitions as a function of magnetic field. 
The normalisation factors (x3 and x30) compensate for a decrease in the transition 
strengths. 

the ground state still not being completely decoupled. For example, the second- 
highest energy state in the ground-level manifold asymptotes to being 1 1/2, 1/2) at 
large field in the \mj,mi) basis. However, at 0.618 T the composition of the state 
is 0.99 1 1/2, 1/2) + 0.14 |— 1/2,3/2). The weak component of this state couples via an 
allowed Amj = —1, Amj = transition to the lowest energy level of the excited 
manifold | — 3/2, 3/2). The theoretical transition strength for the outermost of the weak 
transitions is shown as the solid (yellow) curve for all fields. Numerically diagonalizing 
the atomic Hamiltonian matrix allows one to predict the energy levels and probabilities 
for all the features. The 16 strong transitions (Amj = ± 1) are easily obtained from 
figure [H however, the advantage of our technique is that the detunings and absolute 
linestrengths of the 12 weakly allowed transitions are also given. 
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Detuning (GHz) 

Figure 5. Measured and theoretical absolute dispersion and absorption spectra in 
the hyperfine Paschen-Back (HPB) regime for the Rb D 2 line, through a vapour cell 
(99% 87 Rb, 1% 85 Rb) of length, L = 1 mm, as a function of linear detuning, A/2tt. 
Plot (a) shows solid (red) measured and dashed (black) theoretical transmission in 
the absence of field at a temperature of (127.4 ± 0.8) °C and a width of rbuffer/27r = 
(23.7 ± 1.2) MHz. Plot (b) shows solid (blue) measured and dashed (black) theoretical 
transmission in the presence of a field of (0.618 ± 0.003) T, at a temperature of (114.8 
± 0.6) °C and a width of r huSet /2ir = (47 ± 10) MHz. Below is a plot of residuals 
showing the excellent agreement between theory and experiment in plot (b), with an 
rms of 0.9 %. Plot (c) shows solid (olive) measured and dashed (black) theoretical 
differencing signals in the presence of a field of (0.599 ± 0.003) T at a temperature of 
(126 ± 0.8) °C and a width of r bufier /27r = (63 ± 12) MHz. The bottom plot is the 
residuals for (c), with an rms of 3.3 %. There was no attempt made to add in the shift 
due to the buffer gases. 



5. Dispersion in the hyperfine Paschen-Back regime 

Figure [5] shows the measured and theoretical absolute dispersion and absorption in 
the HPB regime as a function of detuning on the D 2 line in a 87 Rb vapour. Plot (a) 
shows solid (red) measured and dashed (black) theoretical absolute transmission at a 
temperature of (127.4 ± 0.8) °C, showing the Rb absorption features in the absence of 
a magnetic field. Plot (b) shows solid (blue) measured and dashed (black) theoretical 
absolute transmission at a temperature of (114.8 ± 0.6) °C in the presence of a magnetic 
field of (0.618 ± 0.003) T. The transmission signal describing the absorptive properties 
of the medium, was measured using a single photodiode after the cell. Below the 
main plot are the residuals showing excellent agreement over 60 GHz, with an rms of 
0.9 %. Plot (c) shows the solid (olive) measured and dashed (black) theoretical absolute 
dispersion at a temperature of (126.0 ± 0.8) °C in the presence of a magnetic field of 
(0.599 ± 0.003) T. The differencing signal describing the dispersive properties of the 
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medium, was measured using a balanced polarimeter. Below the main plot, residuals 
show good agreement over 60 GHz, with an rms of 3.3 %. In the theoretical and 
measured signals the additional weak features are clearly visible, owing to the fact that 
the ground terms are not completely decoupled. Such a result highlights the strength 
of a model for understanding the energy levels and transition probabilities of such an 
ensemble. In contrast to the zero crossings associated with the dispersive features of 
the off- resonant Faraday effect [9], the various features in figure [5] are associated with 
Zeeman shift resonances and are therefore less sensitive to temperature. Consequently 
this opens up the possibility for locking far off-resonance. 

6. Conclusions 

In summary, we have tested our model for the electric susceptibility of Rb vapour for 
magnetic fields up to and including 0.6 T, which corresponds to the HPB regime on the 
D 2 line. We have demonstrated excellent agreement between the theoretical predictions 
and the experimental measurements of the absolute absorption and dispersion properties 
of the medium. Our study extends the range of magnetic fields for which the 
theoretical model has been tested by an order of magnitude compared to previous 
work. Understanding the optical properties of atomic vapours in the HPB regime 
will find utility in many applications, such as realising compact optical isolators [13J, 
measuring magnetic fields with submicron spatial resolution, and constructing tunable 
atomic frequency references |19j . 
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